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ABSTRACT 

We reassess the estimate of the cross-correlation of the spatial distribution of the NRAO 
VLA Sky Survey (NVSS) radio sources with that of Cosmic Microwave Background (CMB) 
anisotropies from the Wilkinson Microwave Anisotropy Probe (WMAP). This re-analysis 
is motivated by the fact that most previous studies adopted a redshift distribution of NVSS 
sources inconsistent with recent data. We find that the constraints on the bias-weighted red- 
shift distribution, b(z) x Af(z), of NVSS sources, set by the observed angular correlation 
function, w(6), or, equivalently, by the power spectrum of their spatial distribution, strongly 
mitigate the effect of the choice of Af(z). If such constraints are met, even highly discrepant 
redshift distributions yield NVSS-WMAP cross-correlation functions consistent with each 
other within statistical errors. The models favoured by recent data imply a bias factor, b(z), 
decreasing with increasing z, rather than constant, as assumed by most previous analyses. 
As a consequence, the function b(z) x N(z) has more weight at z < 1, i.e. in the redshift 
range yielding the maximum contribution to the IS W in a standard ACDM cosmology. On the 
whole, the NVSS turns out to be better suited for ISW studies than generally believed, even 
in the absence of an observational determination of the redshift distribution. The systematics 
introducing spurious power in the angular correlation function on scales of several degrees 
are strongly reduced restricting the analysis to the sub-sample brighter than 10 mJy. Even 
though this sub-sample comprises less than one third of the NVSS sources, it yields a slightly 
more significant detection of the ISW effect than the full sample (3cr rather than 2.5a), The 
NVSS-WMAP cross-correlation function is found to be fully consistent with the prediction 
of the standard ACDM cosmology. 

Key words: cosmic microwave background — cosmological parameters — cosmology: ob- 
servations — radio continuum: galaxies. 



1 INTRODUCTION 

As first pointed out by Crittenden & Turok (1996) a promising way 
of probing the (linear) Integrated Sachs & Wolfe (1967) (ISW) 
effect is through correlations of Cosmic Microwave Background 
(CMB) maps with tracers of large scale structure. Since the ISW 
effect shows up on large angular scales, a significant application 
of this test had to await until the first all-sky high quality CMB 
temperature maps have been provided by WMAP (Bennett et al. 
2003; Hinshaw et al. 2007). The WMAP data have been cross cor- 
related with a variety of radio, IR, optical, and X-ray surveys (see 
Aghanim et al. 2007 for a review) to look for evidences of a decay 
of the gravitational potential due to the influence of dark energy. 
The first (marginally) significant correlations (at the ~ 2.5a 



level) of the WMAP first year data were found by Boughn & Crit- 
tenden (2004, 2005) with the HEAO-1 A2 full sky hard X-ray map 
(Boldt 1987) and with a map of the number density distribution of 
1.4 GHz radio sources detected by the NRAO VLA Sky Survey 
(NVSS), covering about 80% of the sky. The cross correlation be- 
tween NVSS and WMAP 1st year data was confirmed by Nolta et 
al. (2004), and by Vielva et al. (2006) using different techniques, in 
real, harmonic and wavelet spaces. Further analysis was triggered 
by the release of the WMAP 3 year data. Some groups exploited 
again the NVSS data (Pietrobon et al. 2006; McEwen et al. 2007, 
2008). 

Independently of the analysis technique used, the compari- 
son of the correlations inferred from the data with model predic- 
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tions require two basic ingredients: the redshift distribution and the 
bias parameter of sources. Since redshift measurements are avail- 
able only for a tiny fraction of NVSS sources, it is necessary to 
resort to model estimates. All the investigations quoted above ba- 
sically adopted the redshift distribution given by model RLF1 of 
Dunlop & Peacock (1990) that was shown by Boughn & Crit- 
tenden (2002) to reproduce the autocorrelation function of NVSS 
sources assuming a redshift-independent bias parameter (see also 
Cress & Kamionkowsky 1998, Magliocchetti et al. 1999). How- 
ever, as pointed out by Blake et al. (2004) and Negrello et al. 
(2006), this model overestimates the density of low-redshift radio 
sources (Magliocchetti et al. 2000), which was found to be con- 
sistent with Dunlop & Peacock's (1990) pure luminosity evolution 
model (Magliocchetti et al. 2002). But if the adopted redshift distri- 
bution, and the associated bias model, are incorrect, the conclusions 
of the analysis may be flawed. In fact, the contribution to ISW sig- 
nal in a ACDM cosmology peaks at z ~ 0.4 and rapidly declines 
at higher and lower redshifts (Afshordi 2004). 

This situation has motivated our re-visitation of the problem. 
In § 2 we discuss the interpretation of the autocorrelation func- 
tion of NVSS sources in the light of the results by Negrello et al. 
(2006). In § 3 we provide a new estimate of the WMAP-NVSS 
cross-correlation as a function of the angular separation and of the 
associated errors using both the Internal Linear Combination (ILC) 
temperature map built using the 3 year WMAP data (Hinshaw et al. 
2007) and the foreground cleaned CMB map obtained by Bonaldi 
et al. (2007) using the same data and allowing for the presence 
of an "anomalous" dust-correlated Galactic emission, in addition 
to the canonical synchrotron, free-free and thermal dust emissions. 
In § 4 the observed cross-correlation function (CCF) is compared 
with models prediction. The main conclusions are presented and 
discussed in 5 5. 



2 THE TWO-POINT ANGULAR CORRELATION 
FUNCTION OF NVSS SOURCES 

The NVSS (Condon et al. 1998) is the largest-area radio survey at 
1.4 GHz. It covers ~ 10.3 sr of the sky north of S = -40°; the 
source catalogue contains 1.8 x 10 6 sources. The two-point angular 
correlation function, w(0), of NVSS sources has been measured by 
Blake & Wall (2002a, 2002b) and Overzier et al. (2003) for differ- 
ent flux-density thresholds between 3 mjy and 500 mjy. The overall 
shape of w(6) is well reproduced by a double power-law. On scales 
below ~ 0.1°, the steeper power-law reflects the distribution of the 
resolved components of single giant radio sources. On larger scales 
the shallower power-law describes the correlation between distinct 
radio sources. The latter provides insights on the way in which they 
trace the underlying dark matter distribution, and on the cosmolog- 
ical framework which determines the distribution of dark matter at 
each epoch. In the present context we are interested only on the 
large-scale behaviour. 

Earlier investigations of the ISW effect used NVSS sources 
down to the formal detection limit of ~ 2.5 mjy, although the com- 
pleteness is only ~ 50% at the faintest fluxes, while it rapidly in- 
creases with increasing flux, reaching 99% at 3.4 mjy (Condon 
et al. 1998). Also, systematic surface density gradients, yielding 
spurious contributions to w(9) (or to the power spectrum), are ap- 
proximately negligible only for Si.4GHz *5 10 mjy (Blake & Wall 
2002a). Following Negrello et al. (2006), we will therefore include 
in our analysis only sources brighter than 10 mjy. The NVSS source 
surface density at this threshold is 16.9 deg -2 . 



The angular correlation function of a population of extragalac- 
tic sources is related to their spatial correlation function, £(r, z), 
and to their redshift distribution, Af(z), by Limber's (1953) equa- 
tion: 



w(9) 



dzN 2 (z) \ d(5z)£[r(5z,6),z] 



dzN(i 



(1) 



Here, r (Sz, 6) is the comoving spatial distance between two objects 
located at redshifts z and z + 5z and separated by an angle 9 on the 
sky. 

Equation l[T) shows that the contributions to w(9) of a popula- 
tion of extragalactic sources are weighted by the square of their red- 
shift distribution. The N(z) yielded by Dunlop & Peacock (1990) 
model RLF1, adopted in previous ISW studies, has a spike at low 
z, which, as pointed out by Negrello et al. (2006), is inconsistent 
with the recent accurate determinations of local luminosity func- 
tions of NVSS sources (Magliocchetti et al. 2002; Sadler et al. 
2002; Mauch & Sadler 2007). Also, Ho et al. (2008) found that 
such redshift distribution is inconsistent with the results of cross- 
correlations of NVSS with 2MASS and SDSS samples. As stressed 
by Magliocchetti et al. (1999), because of the Af 2 (z) weighting of 
the projected spatial correlation function, the spike makes a sig- 
nificant contribution to w(8), even though the visibility volume 
of the corresponding sources is small. Actually, such sources, that 
have been interpreted as the spiral and starburst galaxy populations 
emerging at mJy levels, would provide the dominant contribution 
to the clustering signal on angular scales larger than ~ 2° (Cress & 
Kamionkowski 1998; Magliocchetti et al. 1999; Blake et al. 2004). 
But the contribution of spiral and starburst galaxies to counts at 
5i.4GHz ^ lOmJy, computed using the recent determinations of 
their luminosity function, is very small (~ 0.5%; Negrello et al. 
2006). Yet, a statistically significant positive w(8) was measured 
on scales ^ 5°, and the interpretation of this data proved to be 
challenging. 

On scales large enough for the clustering signal be due to 
galaxies residing in distinct dark matter haloes and under the as- 
sumption of a one-to-one correspondence between sources and 
their host haloes, the spatial two-point correlation function can be 
written as the product of the correlation function of dark matter, 
(dm, times the square of the bias parameter, b (Matarrese et al. 
1997; Moscardini et al. 1998): 



£(r,z) = b (Af ff , z)£dm(»", z) 



(2) 



where M c a represents the effective mass of the dark matter haloes 
in which the sources reside. The observationally determined spa- 
tial correlation length of star-forming galaxies, ro(z — 0) ~ 3— 
4Mpch _1 (Saunders et al. 1992; Wilman et al. 2003) is consistent 
with an effective halo mass not exceeding Mjg = 10 11 Mq h _1 . 
Together with the small contribution of these sources to the counts 
for Si.4GHz lOmJy, this implies that they are negligible con- 
tributors to the observed w{0) (see Negrello et al. 2006 for more 
details). 

The clustering properties of low-z AGN-fuelled radio sources 
are much stronger and imply that they reside in dark matter haloes 
more massive than ~ 10 13 - 4 M Q (Magliocchetti et al. 2004). 
These could then easily account for the observed w(9) on scales 
<; i°_2° . However, a generic prediction of the Cold Dark Matter 
paradigm of structure formation, is that the spatial correlation func- 
tion of matter displays a sharp cut-off around a comoving radius of 
a few hundred Mpc, and becomes negative on larger scales. For 
typical redshifts z ~ 1, this translates in a cut off of w{6) on scales 
~ 2° . On larger scales, the negative contribution of z ~ 1 sources 
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overwhelms the positive contributions from low-z sources if M c r 
is redshift-independent, as indicated by data on optically selected 
quasars (Porciani et al. 2004; Croom et al. 2005). To account for 
the data in the framework of the standard hierarchical clustering 
scenario it is necessary that the effective bias factor decreases with 
increasing z, rather than strongly increase like for optical quasars. 
This adds weight to the contribution of low-z NVSS sources to the 
matter density fluctuations, enhancing the expected CMB-NVSS 
cross-correlation. 

Negrello et al. (2006) found that the observed w(8) can be 
reproduced if M c g is proportional to the characteristic mass of 
virialized systems, M*(z), with M cfi (z = 0) ~ 10 15 M . This 
implies M c ff (z = 1.5) ~ 10 13 Mq, close to the value found for 
optical quasars. A possible interpretation of this result is that AGN- 
powered radio sources are so rare locally that the typical halo host- 
ing one of them must include a very large number of galaxies, i.e. 
must correspond to a rich cluster of galaxies. The strong evolution 
of radio sources implies that they become relatively more abundant 
with increasing redshift, and the typical halo mass decreases corre- 
spondingly. This happens in such a way that the clustering of radio 
sources reflects that of the largest haloes that collapse at any cosmic 
epoch. 

We have repeated the analysis using the 2 versions (MEAN-z 
and HIGH-z) of the luminosity/density evolution (LDE) model by 
Dunlop & Peacock ( 1 990), which are also consistent with the recent 
data on the local luminosity function. As illustrated by Fig. [TJ the 
requirement that the observed w(8) is reproduced (lower left-hand 
panel) forces (lower right-hand panel) the bias-weighted redshift 
distributions, b(z) x N(z), corresponding to the 3 models to be 
very close to each other in the redshift interval from which most of 
the contribution to the ISW signal comes, thus lessening the effect 
of uncertainties in the redshift distribution and on the bias factor 
(upper panels). Also the redshift dependence of the bias factor adds 
weight to the interesting redshift range, implying that the NVSS 
sample is very well suited for ISW studies. 



3 OBSERVATIONAL ESTIMATE OF THE WMAP-NVSS 
CROSS-CORRELATION FUNCTION 

The empirical cross-correlation function was computed as: 

(fc) ~ E« ' 

where Ni is the number of NVSS sources in the i-th pixel of the 
map of number count fluctuations, Tj is the temperature fluctua- 
tion in the j-th pixel of the CMB map, the two pixels being at an 
angular distance 8 k - The weights wj and wf 1 are equal to unity 
for valid pixels or equal to zero if the pixels fall within a region 
not covered by the survey or masked; the apex T and N refer to 
the CMB and to the NVSS maps, respectively. The maps are pix- 
elized with the Healpix (Gorski et al. 2005) resolution parameter 
NSIDE=32, corresponding to a pixel size of about 1°.8 x 1°.8. 

We used two CMB maps: the WMAP three-year Internal 
Linear Combination (ILC) map (Hinshaw et al. 2007) and the 
foreground-cleaned map yielded by model M2 of Bonaldi et al. 
(2007). As for the NVSS map, our choice of a 10 mJy flux limit has 
dispensed us with the need of introducing corrections for the spuri- 
ous density gradients in the survey noted by Boughn & Crittenden 
(2002) for the sample limited to 2.5 mJy. We have removed the 
strip at \b\ ^ 5° where extragalactic sources are blurred by bright 
Galactic sources and by small scale structure of the synchrotron and 



free-free emission. Including the unobserved region at S < —40°, 
the NVSS mask covers 25% of the sky. 

The CMB mask adopted is based on the KpO mask (Bennett et 
al. 2003), which excludes the regions most contaminated by Galac- 
tic foregrounds and the brightest point sources in the WMAP bands. 
The KpO mask has NSIDE=512, corresponding to a pixel size of 
about 7'. We degraded this map to the lower resolution adopted for 
this work and we set wf — for all pixels including at least one 
masked 7' pixel. The resulting mask excludes 37% of the sky. The 
union of NVSS and CMB masks covers 48% of the sky. 

For both CMB maps, we computed the empirical CCF 
[eq. l[3}] for a set of angular separations, in degrees: 8k = 
{0.0, 1.8, 3.6, 5.4, 7.2, 9.0, 10.8}. As discussed by Blake & Wall 
(2002a) the tiny angular correlations among NVSS sources on 
scales > 5° [w(8) ~ 10~ 4 ] may well be due to small system- 
atic variations in the source surface density due to calibration prob- 
lems. This may also add some spurious power to the CCF; the 
effect is however small and the empirical CCF is indeed found 
to be consistent with zero on these scales. For each 8k we con- 
sidered pixel pairs whose centers have separations in the range 
8 k -A8<8<8 k + A8 with A8 = 1°. To estimate the errors we 
simulated 1000 mock NVSS maps by randomly redistributing the 
unmasked pixels of the true NVSS map. For each mock NVSS map 
we computed the CCF with the CMB map. The error on CCF(#k) 
is then estimated as the rms value of the CCFs for the mock maps. 
Alternatively, the significance of the signal can be estimated from 
the fraction of mock CCFs at or above its amplitude. We find that 
the significance estimated in this way is in excellent agreement with 
that expected from the S/N ratio for a Gaussian distribution. The re- 
sults are shown in Fig. [2] The significance of the measured signal 
is similar for the two CMB maps used. The CCF differs from zero 
at the 2.1-2.4(7 level for each of the first 2 points, and the overall 
significance of the ISW signal is 99.7% or slightly higher, corre- 
sponding to 3a or somewhat more. Boughn & Crittenden (2003) 
reported a somewhat less significant (2.5<r) detection of the NVSS- 
CMB cross-correlation, suggesting that the systematics affecting 
the distribution of the faintest NVSS sources overcome the statisti- 
cal advantage of their far larger number. 



4 THE THEORETICAL CROSS-CORRELATION 
FUNCTION 

Following Nolta et al. (2004) we write the cross-correlation power 
spectrum between the surface density fluctuations of NVSS sources 
and CMB fluctuations as: 

Cr = (aLaL) = 4^"^ fA\k)f?(k)tf(k), (4) 

where /j and ff are the NVSS and CMB filter functions, respec- 
tively, and A 2 (fc) is the logarithmic matter power spectrum today 
as a function of the wave number k: 

1.3 / rk \ 3+n 

= ^ p ^ k ) = (£) r /W- © 

For Sh we have used the expression given by eq. (A3) of Eisenstein 
& Hu (1998), valid for a flat universe with a cosmological con- 
stant. The reference values of cosmological parameters are those 
given in Table 2 of Spergel et al. (2007), namely Ho = 100ft = 
72kms~ 1 Mpc~ 1 ,f2 m = 0.13/i~ 2 ,ft A = l-d m ,n = 0.96, and 
erg = 0.78. The constant A was computed normalizing the power 
spectrum to as = 0.78; we find A = 1.06. The matter transfer 
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Figure 1. Redshift distributions (upper left-hand panel), evolution of the bias factor (upper right-hand panel), bias-weighted redshift distribution (lower left- 
hand panel), and fits of the observed angular correlation function (lower right-hand panel) of NVSS sources brighter than 10 mJy for Dunlop & Peacock (1990) 
pure luminosity evolution (PLE) model (solid line) and the 2 versions of luminosity/density evolution (LDE) model, MEAN- 2 (dashed line) and HIGH-2 (dot- 
dashed line). The three model w{9) include the small systematic offset &w(6) = 10 -4 attributed to small systematic variations in the source surface density 
due to calibration problems (Blake & Wall 2002a). The three dots-dashed lines show, for comparison, the results obtained with Dunlop & Peacock (1990) 
RFL1 model with a constant value of the bias parameter (b = 1.7), as assumed by most previous analyses. The data points with error bars show the Blake & 
Wall (2002b) estimates of w{6). 



function, Tf(k), was computed using CMBfast (Seljak & Zaldar- 
riaga 1996). 

The NVSS filter function writes: 



f?(k) 



dz ^ bAz) m jl{kcil{z)) > 



(6) 



where (dN/dz)dz is the mean number of sources per steradian 
with redshift z within dz, brighter than the flux limit, b r (z) is the 
bias factor relating the source overdensity to the mass overdensity, 
assumed to be scale-independent, D(z) is the linear growth fac- 
tor of mass fluctuations, ji (x) is the spherical Bessel function, and 
rj(z) is the conformal look-back time: 



7/0) = 



dz' 
H(z) 



dz' 



HoE(z') 



(7) 



with E(z) = y/^A + + z) 2 + n m (l + z) 3 , n R being the 

curvature parameter, set to zero in this paper. For the linear growth 
factor, appropriate for the large scales of interest here, we adopt the 
approximation (Carroll et al. 1992): 



D{z) = WT^){ nm{z)4/7 ' nA{z) + 



1 + 



n A {z) 



70 



(8) 



withQ m 0) = n m (l + zf/E(z) 2 and Q A (z) =ft A /E{z) 2 . 

As mentioned in § 1, we restrict ourselves to NVSS sources 
brighter than 10 mJy and adopt the redshift distributions given by 
three Dunlop & Peacock (1990) evolution models for AGN pow- 
ered sources. As for star-forming galaxies, we estimate the redshift 
distribution using the local luminosity function of Magliocchetti et 
al. (2002) and no evolution: evolutionary effects, if any, are irrele- 
vant for our sample since star-forming galaxies brighter of 10 mJy 
are found to be at z ^ 0.1. The filter function for the ISW effect is: 



d[D(z)(l + z)/D(0)] 



dz- 



dz 



(9) 



Since eq. ((9) holds for linear perturbations, i.e. for k <C 1, the in- 
tegration over k of eq. ® requires some caution. We have checked 
for a set of multipoles I that the cross-correlation function ob- 
tained using the power spectrum of CMB temperature perturba- 
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Figure 2. Estimates of the WMAP-NVSS cross-correlation function (CCF) 
using the ILC map (triangles) and the foreground cleaned map by Bonaldi 
et al. (2007; diamonds), with lcr error bars. The ILC results are slightly 
shifted to the right to improve the readability. The empirical CCF is com- 
pared with expectations from the best-fit WMAP 3-year cosmology for the 
models described in the text. As in Fig.[T] the solid line refers to the PLE 
model, the dashed line to the LDE MEAN-2 model, the dot-dashed line to 
the LDE HIGH-2 model, and the three dots-dashed line to the model used 
by Nolta et al. (2004). 



tion given by CMBfast are well approximated by setting in eq. l[4} 
fc m ax = 0.02 Mpc" 1 . 

The cross-correlation function as a function of the angular sep- 
aration 8 is then obtained as: 



c NT (e) = y 2 i±± c r Pilose). 



(10) 



where P\ are the Legendre polynomials. For the range of scales of 
interest here, it is enough to sum over the interval 2 ^ I ^ 200. 

Figure [2] shows that, once the constraints from the observed 
w(0) are taken into account, the predicted CCF for the WMAP 3- 
year values of cosmological parameters is only weakly affected by 
the uncertainties on the redshift distributions of NVSS sources. We 
find a good agreement between the predicted and empirical CCFs, 
somewhat at odds with the finding by Ho et al. (2008) that the 
WMAP 3-year model predicts an ISW amplitude about 2a below 
their estimate, which exploits various galaxy surveys in addition 
to the NVSS. The strong constraints from the angular correlation 
function of the latter sources are not taken into account, however. 



5 DISCUSSION AND CONCLUSIONS 

The NVSS sample is of primary importance for the investigation of 
the ISW effect because of its very large area (82% of the sky), its 
large number of sources (almost 2 x 10 6 , its depth in redshift, and 
its immunity to systematics related to the uncertain corrections for 
dust extinction. It has however drawbacks, not all of which were 
properly dealt with by previous analyses. On one side, spurious 
contributions to the power spectrum of the NVSS source distribu- 
tion is introduced by known systematic effects such as small sys- 
tematic variations in the source surface density due to calibration 
problems at low flux densities and to striping effects, particularly at 
negative declinations (Blake & Wall 2002a; Boughn & Crittenden 
2002). To minimize this problem, we have restricted our analysis 



to sources with 51.4GHz ^ 10 mJy (Blake & Wall 2002a,b); at this 
flux limit, striping effects are negligible. 

The main potential problem, however, stems from the very 
limited direct redshift information, forcing us to resort to models. 
Most previous analyses (Boughn & Crittenden 2002, 2004; Nolta 
et al. 2004; Pietrobon et al. 2006; Vielva et al. 2006; McEwen et al. 
2007) relied on model RLF1 of Dunlop & Peacock (1990), which 
was found by Boughn & Crittenden (2002) to reproduce the NVSS 
autocorrelation function. However, as pointed out by Negrello et al. 
(2006), this model is inconsistent with the recent determinations of 
the luminosity function of low-z NVSS sources (Magliocchetti et 
al. 2002; Sadler et al. 2002; Mauch & Sadler 2007). A totally inde- 
pendent confirmation of this conclusion has been obtained by Ho et 
al. (2008) by cross-correlating the NVSS with 2MASS and SDSS 
samples, whose redshift distributions are known. In this way these 
authors obtained constraints on the bias-weighted redshift distribu- 
tion, b(z) x N(z), which is the relevant quantity for predicting the 
NVSS-CMB cross-correlation induced by the ISW effect [eq. 

The large spread of redshift-dependent radio luminosity func- 
tions yielded by the original set of Dunlop & Peacock (1990) mod- 
els is substantially narrowed down by the more recent data on the 
radio luminosity function at low redshifts, that are fully consis- 
tent with the pure luminosity evolution and by the two luminos- 
ity/density evolution models. Strong constraints on b(z) x N(z) 
follow from the observational determination of the angular correla- 
tion function, w(0), of NVSS sources, which depends on [b(z) x 
Af(z)] 2 . In fact, the functions b(z) xAf(z) for three models consid- 
ered here are found to be substantially closer to each other, particu- 
larly at low-redshifts, than the corresponding redshift distributions. 
The decrease of the bias factor with increasing z implied by obser- 
vational determinations of w(6) adds weight to the redshift range 
yielding most of the contribution to the ISW signal, making the 
NVSS sample even better suited to test the effects of dark energy 
on the growth of structure than one would guess from the estimated 
redshift distribution. 

A concern expressed in some previous papers is that the mi- 
crowave emission of radio sources themselves contribute to the ob- 
served WMAP-NVSS correlations. However, the radio source con- 
tribution to temperature fluctuations is dominated by the brightest 
sources, while the contribution to the source number density fluc- 
tuations, relevant for studies of the ISW effect is dominated by 
the faintest sources. Number density and surface brightness fluc- 
tuations are thus largely uncorrelated, implying that the source mi- 
crowave emission rather weakens the correlations. We find that, if 
no sources are masked, the CCF estimates decrease by 6-10%. On 
the other hand, removing just the WMAP sources (Hinshaw et al. 
2007; Lopez-Caniego et al. 2007) is enough to make the contami- 
nation of the ISW signal negligibly small, consistent with previous 
results (e.g. McEwen et al. 2007; Ho et al. 2008). 

In summary, the main results of this paper are: 



• the constraints on the bias-weighted redshift distribution, 
b(z) x Af(z), of NVSS sources, set by the observed angular cor- 
relation function or, equivalently, by the power spectrum of their 
spatial distribution, strongly mitigate the effect of the large uncer- 
tainties on N(z). Even highly discrepant redshift distributions, in- 
cluding instances ruled out by recent data, yield NVSS-WMAP 
cross-correlation functions consistent with each other within statis- 
tical errors, although the use of the Dunlop & Peacock (1990) RLF1 
model leads to a slight underestimate of the predicted ISW signal 
for the WMAP 3 year values of the cosmological parameters. This 
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fact largely alleviates what was perceived as the main weakness of 
the exploitation of the NVSS for studies of the ISW. 

• The models favoured by recent data imply bias functions, 
b(z), decreasing with increasing z, rather than constant, as as- 
sumed by most previous analyses. As a consequence, the func- 
tion b(z) x N(z) has more weight at z < 1, i.e. in the redshift 
range yielding the maximum contribution to the ISW in a stan- 
dard ACDM cosmology. Again, this means that the NVSS is better 
suited for ISW studies than generally believed. 

• The systematics (striping, incompleteness, etc.) introducing 
spurious power in the angular correlation function on scales of sev- 
eral degrees are strongly reduced restricting the analysis to the sub- 
sample with a flux limit of 10 mjy. Even though this sub-sample 
comprises less than one third of the NVSS sources, it yields a 
slightly more significant detection of the ISW effect (3<j rather than 
2.5a). The microwave emission of NVSS sources weakens the cor- 
relations with WMAP fluctuations, rather than increasing them as 
previously argued. The effect is, however, negligibly small. 

• The NVSS-WMAP cross-correlation function is found to be 
fully consistent with the prediction of the standard ACDM cosmol- 
ogy- 
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